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Realistic selfconsistent calculations of unconventional superconductivity in a heavy-fermion mate- 
rial are reported. Our calculations for UPdgAls start from accurate energy band dispersions that 
are computed within the local spin-density functional theory and provide Fermi surfaces in agree- 
ment with experiment. Using physically motivated, realistic pairing potentials it is shown that the 
superconducting gap has two lines of nodes around the z-axis, thus exhibiting d-wave symmetry in 
the Aig representation of the Dgh point group. Our results suggest that in a superconductor with 
gap nodes, the prevailing gap symmetry is dictated by the constraint that nodes must be as far as 
possible from high-density areas. 



PACS numbers; 74.25. Jb, 74.20.-z 

Superconductivity (SC) in heavy-fermion (HF) mate- 
rials exhibits a fascinating complexity of phenomena [Q . 
The nature of the SC state is mostly anticipated to be 
unconventional, i.e., an additional symmetry is broken 
in the SC state, and the order parameter is therefore 
not the conventional spin singlet s-wave type (2|j^. The 
identification of the order parameter symmetry is an es- 
sential issue in current studies of unconventional SC. A 
common approach to address the symmetry of the or- 
der parameter is through investigations of the SC gap, 
which is expected to have the same symmetry as the or- 
der parameter. While progress has been made for high- 
Tc superconductors (HTSC), the identification of the gap 
symmetry remains an open question for most HF materi- 
als. Unfortunately, the eventual identification of the gap 
symmetry in a HF compound is not sufficient to identify 
the pairing mechanism. Additional elements are neces- 
sary and in the case of HF materials, it was suggested 
that spin-fluctuations - either in the itinerant or lo- 
calized limit as magnetic excitons |^,^ - may replace the 
phonons as mediators of the pairing. 

In order to identify the gap symmetry and clarify its 
relationship to the pairing mechanism, selfconsistent so- 
lutions of the SC gap equations with physical momen- 
tum dependent pairing potentials and accurate band 
structures are unavoidable. Such type of calculations 
have been achieved only recently for HTSC and ruthen- 
ates where usually two-dimensionality is assumed and 
a tight-binding dispersion of one or two electron bands 
is considered ||^-|Tl|]. In the case of HF's the situa- 
tion is much more complex. In addition to the three- 
dimensionality, we must deal with several anisotropic 
bands that cross the Fermi level {Ep) producing numer- 
ous highly anisotropic Fermi surface (FS) sheets. Due to 
the complicated structure of the bands in HF materials, 
simple tight-binding fits are impossible to obtain. As a 
consequence, no realistic selfconsistent calculations of SC 



have been achieved so far in HF materials. 

Here we report for the first time a direct combination 
of relativistic band-structure calculations with realistic 
selfconsistent calculations of SC in a HF material. We 
focus on UPd2Al3 which is a fascinating HF supercon- 
ductor having a moderately large specific heat coefficient 
7 = 140mJ/molK^, and a relatively high critical temper- 
ature Tc = 2 K 1 12 1 . It orders antiferromagnetically below 
Tjv — 14.3 K with an ordered moment of 0.85 /xs/U-atom 
p3[ , which is large compared to the moments of other HF 
superconductors as, e.g., UPta, which has a moment of 
only 0.03 /is . A further anomalous feature is that the 
antiferromagnetic (AFM) order coexists with SC below 
2K. Coexistence of magnetism and SC was observed for 
other materials, e.g., containing 4/ elements, but in those 
cases the magnetism is due to the localized 4/ electrons 
far from Ep, whereas the SC is carried by itinerant elec- 
trons at Ep. Conversely, for UPd2Al3 most of the recent 
studies reveal that the SC, magnetic order, as well as 
HF behavior all involve the uranium 5/ states |l^ , |l5| . 
Moreover, some striking similarities with cuprate HTSC 
emerged recently. Neutron scattering experiments [|l6[|l7[ 
have reported a distinct resonant mode that develops be- 
low Tc as in HTSC. Tunneling spectroscopy experiments 
|l8| in the SC state reported peculiar dip-hump features 
in the tunneling spectrum, similar to those observed in 
HTSC. 

We solve selfconsistently the simplified gap equation 
starting from energy dispersions which we computed 
within the framework of density-functional theory in the 
local spin-density approximation (LSDA) . In our compu- 
tational approach we adopt several assumptions. First, 
we approximate the multiband gap equation [l9| by an 
effective singleband equation. This is justified because 
the bands at Ep have mainly 5/ character [Q, and 
there is thus no need to discriminate different band char- 
acters. The pairing potential is therefore independent 
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of the band character. In the multiband Hamiltonian 
pof the fermion operators for the different band charac- 
ters become thereby identical, and an effective singleband 
Hamihonian, with the energy dispersion ek replaced by 
X^n ^nk (n is the band index), results. We obtain con- 
sequently one gap symmetry defined for the whole FS. 
This corresponds to the experimental studies that iden- 
tify a global gap symmetry and not one for every FS sheet 
| p5| , p^j20| . Second, since it is experimentally firmly es- 
tablished that UPd2Al3 is a spin singlet superconductor 
pl| , we explore in our calculations only singlet states. 
Third, we adopt an adiabatic approach in our investiga- 
tion of the SC order parameter. This means that we ne- 
glect the renormalization of the AFM bands by the SC, 
which is justified because the magnitude of the energy 
scale of SC is an order of magnitude smaller than that of 
the AFM order. Fourth, in UPd2Al3 there exist certainly 
many-particle effects that are esponsible for the high ef- 
fective mass, which can cause a renormalization of the 
energy bands. Such processes are not taken into account 
in our LSD A calculations. Previous studies demonstrated 
a very close correspondence of the measured and cal- 
culated de Haas- van Alphen frequencies, and conse- 
quently the FS of UPd2 AI3 is accurately described by the 
LSD A energy dispersions. Since the relevant physics for 
SC happens in the vicinity of the FS, the LSDA bands 
are the appropriate starting point for the study of the 
gap symmetry. Moreover, we focus here on a qualita- 
tive analysis of the gap symmetry avoiding quantitative 
aspects of SC for which the incorrect values of the effec- 
tive masses in our LSDA scheme will indeed have some 
influence. 

We computed the energy bands of UPd2Al3 with the 
relativistic augmented-spherical-wave method p^ . The 
resulting LSDA band structure of UPd2 AI3 along charac- 
teristic symmetry lines is shown in Fig. 1 for the param- 
agnetic and AFM state. The FS resulting from our AFM 
bands consists of four types of sheets, and is practically 
identical to that computed previously jlj] , as is also the 
ordered AFM moment of 0.81 /is, which agrees with the 
experimental moment of 0.85 /ib fl^ . 

The simplified gap equation in the T = regime is 

Ak = i5]y(k,k')AkKXk- + A^,)-l/^ (1) 

k' 

where l^(k, k') is the pairing potential, and A"k' = 
(Sn Cnk)^ is the supcrposition of all the relevant bands 
(dotted lines in Fig. 1). The momentum summation is 
over the full Brillouin zone (BZ). An unphysical property 
of Eq. (1) may occur when summed energies are zero, but 
not the individual energies. To avoid such situation wc 
approximate « e^^.. We verified that this approx- 
imation did not affect the resulting gap symmetry. We 
consider two types of pairing kernels that were shown to 
produce unconventional order parameters. The small-g 



pairing potential adopts at small wavevectors the follow- 
ing form in momentum space l^(k, k') = ^3q:^E^pp"+ 
/i*(k — k'). This kernel is characterized by a smooth 
momentum cut-off Qc which selects the small-q processes 
in the attractive part |p|,p^-p7[ . At larger wavevectors 
the repulsive Coulomb pseudopotential /i*(k — k') pre- 
vails. This type of attractive interaction may be due 
to phonons if screening with short range Hubbard-like 
Coulomb terms is involved ||2^ . It has been claimed that 
this type of kernel may account for unconventional SC 
in HTSC ||,|2^ and in other materials, including HP's 
p5| . Agterberg et al. |^ have recently suggested for 
HF superconductors a multipocket FS with an attractive 
intrapocket and repulsive interpocket potential. Such sit- 
uation results naturally from a potential dominated by 
small-q attractive pairing as the one considered here, if 
the FS contains multipockets. We also considered the 
case of pairing mediated by spin-fluctuations using a 
phenomenological Millis-Monien-Pines pairing potential 
p9| which has the form V^(k, k') = q:j_(.(k-^'-~Q)ii '^liere 
Q = (0,0, 7r/2c) in UPd2Al3 (c is the c-axis lattice con- 
stant). This type of pairing has been suggested for virtu- 
ally all unconventional SC's, including UPd2Al3 | |3C| , p^ . 

To obtain selfconsistent gap solutions with such mo- 
mentum dependent kernels we have used a Fast-Fourier- 
Transform technique. The problem is solved iteratively 
on a symmetric part of the BZ that we discretize with a 
128 X 128 X 128 momentum grid. For each point of our 
3-D grid wc have our LSDA energy bands as an input. 
Within our procedure the momentum space problem is 
fully resolved numerically without any simplification or 
bias on the resulting gap symmetry. 

In UPd2Al3, NMR measurements revealed a singlet SC 
state with a gap of even parity. Because UPd2Al3 
and thus our LSDA bands obey the hexagonal Dq^ point 
group symmetry, the only even parity accessible gap 
states with nodes should transform according to the irre- 
ducible representations of D^h [0 shown schematically 
in Fig. 2. As a first step in our iteration scheme, the 
initial gap configuration is chosen randomly and the sys- 
tem is totally free to converge within the iteration cycle 
towards the most favorable representation. In a second 
step we bias the system by adopting as initial gap con- 
figuration each of these representations and try to force 
convergence towards the chosen representation. In that 
way, we check the ground state character of the solution 
obtained in the first step. The precision of our numeri- 
cal calculation is such that it allows to distinguish which 
representation corresponds to the lowest free energy. The 
lowest-in-energy representation, e.g., has a free energy 
that is 10% less than the next higher one. The second 
step in addition allows us to identify any eventual degen- 
eracy between the various accessible gap symmetries. 

The central result of our calculations is that the only 
accessible gap states with nodes in UPd2Al3 transform 
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according to the fully symmetric irreducible representa- 
tion Aig. This is true irrespective of the precise momen- 
tum structure of the pairing kernel. We show in Fig. 3 
some examples of our selfconsistently obtained gap func- 
tions in the plane L — A — L — L — M — L which contains 
the z-axis (the A — F — A axis) and in the plane ob- 
tained by a 7r/6 or 7r/2 rotation around the z-axis (the 
plane H — A — H — H — K — H). All selfconsistent solu- 
tions shown in fig. 3 have two lines of nodes practically 
perpendicular to the A — F — A axis, thus belonging in 
the Aig representation. No other representation possess- 
ing nodes was accessible for all parameters in the pairing 
kernels that we have investigated. To study the influ- 
ence of the coexisting AFM order on SC we have made 
calculations using both paramagnetic and AFM bands. 
In both paramagnetic and AFM cases, the only solutions 
with gap nodes are of Aig type. Note that our results 
with the spin-fluctuations kernel confirm the model anal- 
ysis by Huth et al. |§. 

There is overwhelming experimental evidence that in- 
deed UPd2Al3 has a gap with the Aig node structure. 
Recent tunnel measurements along the z-axis showed the 
absence of a node in this direction However, the 

presence of nodes is definitely established in UPd2Al3, 
and since Aig is the only even parity gap with nodes 
which is nodeless along the z-direction, it was deduced 
a posteriori that the gap symmetry must be Aig type 
p8[ . In addition, the observation of a spin-fluctuations 
peak below Tc at (0, 0,7r/2c) 17| would be forbidden 
by the involved coherence factors unless the gap changes 
sign along the z-axis as in the Aig representation [ p2[ . 
Finally, measurements of the angular dependence of the 
critical field indicate two lines of nodes perpendicular 
to the z-axis |Q again suggesting the Aig representa- 
tion. UPd2Al3 is perhaps the only HF superconductor 
for which such a clear picture of the nodal structure of the 
gap is obtained from the experiments. We consider the 
agreement of our results with the experiments as strong 
support of our calculations. 

The surprising robustness of the Aig solution results 
from the following general rule: If a system must choose 
a gap symmetry with nodes because of the repulsive ef- 
fective interaction at large wavevectors (short distances) 
it chooses the representation that has the minimum num- 
ber of nodes as far as possible from the high-density ar- 
eas. In fact, the system must maximize the condensa- 
tion free energy and this is obtained when there is a gap 
in the high-density areas and the node (gapless) areas 
are kept minimal. Examining the paramagnetic bands 
of UPd2Al3 (Fig. lb) we see that the high-density areas 
near the FS due to saddle points in the band dispersions 
are found essentially near the A point and near the H 
point. The Aig representation prevails because it is the 
only one without a node near the A point. In the case of 
the AFM bands (Fig. Ic), the saddle points of the bands 
near the FS are found essentially in the z = plane (near 



the r point, along the F — K, K — M and M — F symme- 
try lines) and near the F point along F — A' (A' = A/2). 
The high-density areas in the z — plane exclude the 
Big , Eig and B2g representations which have a line node 
in the z — plane. From the remaining Aig, E2g and 
A2g representations, Aig is the only one without a line of 
nodes parallel to the A — F — A axis which would cross the 
high-density areas at F and near F along F — A'. Also, 
due to the AFM symmetry a high density occurs again 
at A, as A is equivalent to F. The node perpendicular to 
the F — A path in the Aig representation is close to A' 
and therefore does not cross the high-density areas near F 
and A. Note that A2g: Big and B2g are also handicapped 
by the fact that they have more nodal areas. 

To stress the generality of this argument we reconsid- 
ered the case of cuprate HTSC. Saddle points produce 
high-density areas in the vicinity of (0, tt/o) and symme- 
try related areas (a is the lattice constant in the basal 
plane). The dx2-y2 gap symmetry corresponds to nodes 
along the (tt/o, tt/o) direction as far as possible from the 
high-density (0, ±7r/a) and (±7r/a,0) areas. Using both 
types of pairing kernels as adopted for UPd2Al3 (now 
Q = (tt/c, 7r/a)) we were never able to obtain a dxy solu- 
tion for which the nodes cross the high-density (0, 7r/a) 
areas no matter the details of the interaction. 

In conclusion, we have computed the first selfconsis- 
tent solutions of SC in the HF material UPd2Al3. Our 
calculations demonstrate that the only even parity ac- 
cessible gap symmetry with nodes transforms according 
to the Aig irreducible representation of the Dg/i point 
group, independent of whether the pairing potential is 
dominated by small-q attractive processes or by spin- 
fluctuations. The robustness of the Aig representation 
is analyzed to be due to the presence of high-density re- 
gions of the relevant bands in the vicinity of the F and 
A points, that preclude node formation at these points. 
As a general rule we obtain that nodes must be as far as 
possible from high-density areas in the phase space. 
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F. Stcglich for illuminating discussions. We thank M. 
Huth for the figure of the irreducible representations of 
the Deh group. 
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FIG. 1. (a) The hexagonal Brillouin zone with high sym- 
metry points. In the AFM phase the BZ is reduced by a factor 
two along the z-axis (denoted by the primed letters), (b) The 
energy bands along high symmetry directions in the param- 
agnetic, and (c) the AFM phase of UPd2Al3. The relevant 
bands for superconductivity are shown by the dotted lines. 



FIG. 2. The even parity irreducible representations of 
the Dsh point group. While in principle all are accessible 
for the gap, only the A\g results from our selfconsistent gap 
calculations. 



FIG. 3. (color) Some examples of the selfconsistent gap so- 
lutions in characteristic planes (with sides indicated) that con- 
tain the r point (in the center) and the z-axis (A — F — A axis). 
The thick black lines show the nodes and in the blue areas the 
gap is negative. All nodes cut the A — F — A axis as in the A\g 
representation. We show here results obtained with small-q 
pairing and AFM energy bands, (a and b), spin-fluctuations 
pairing and AFM bands (c), and small-q pairing and param- 
agnetic bands (d). In (d) the sign of the gap is inverted to 
visualize better the variation of the gap. 
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